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ABSTRACT 

A thin, flat, elastic plate of length L and infinite 
width lies in the x-y plane. The rest of the x-y plane is 
occupied by a rigid baffle that separates the upper z space 
from the lower. A harmonic force drives the plate from 
below, and a compressible air stream with a viscous boundary 
layer flows parallel to the upper surface along the length L. 
Linear partial differential equations governing the forced 
response of the coupled plate-aerodynamic system are derived 
along with appropriate boundary conditions. Solution of 
these equations involves modal series solution of the plate 
equation, Fourier transform of the aerodynamic equation in 
the direction parallel to the plate, and Frobenius - series 
or Runge-Kutta numerical solution of the resulting o.d.e. 
with variable coefficients (due to the boundary layer) in 
the direction normal to the plate. The aerodynamic solution 
depends on the dimensionless ratios Mach number, Prandtl 
number, temperature boundary condition, boundary layer 
thickness and acoustic coincidence number. The plate 
response solution depends in addition on the parameters 



plate frequency ratio, density coupling number, plate 
eigenvalue spectrum, and the force distribution. Numerical 
results for Mach number, M, of 0.5 are presented to 
illustrate the method of solution and the effects of boun- 
dary layer thickness on the fluid pressures along the plate 
and on the plate response. Calculations for limiting cases 
are presented showing that the method produces appropriate 
results for M = 0 , and for M = 0.5 but zero boundary 
layer thickness. The surface pressure distributions show 
that peak pressure shifts forward and decreases slightly as 
the boundary layer thickness increases. Variations of 
boundary layer thickness caused large percentage changes of 
the "virtual mass" and aerodynamic damping associated by 
analogy with the fluid effects on the plate dynamic response. 
Calculations of basic solution parameters for a linear 
velocity profile and for a Blasius profile showed that the 
same system response could be obtained from each profile 
provided that an appropriate value of boundary layer thick- 
ness was chosen for each profile. 
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Chapter I 
Introduction 


The practical problem 

Turbulent boundary layers are often formed on the sur- 
faces of aircraft in flight. The response of the aircraft skin 
structure to the fluctuating pressures associated with the 
turbulence can lead to structural fatigue failure due to 
fluctuating stress response, and can lead to unwanted 

transmission of noise from the boundary layer to the interi- 

(2 2 ) 

or of the aircraft. ' The ability to predict the struc- 
tural responses is required in order to design structures 
that will withstand fatigue and that will minimize the trans- 
mission of noise. 

Panel- Aerodynamic Coupling 

The response of aircraft structural panels is affected 
in an important way by its coupling to the passing air 
stream, in addition to the fact that the turbulence in the 
air stream is driving the panel motions. When the panel 
responds to the turbulence, its motion generates additional 
pressures in the air stream; these pressures act on the panel 
to alter its motion. The panel motion is thus coupled to the 
pressure in the air stream in such a way that calculation of 
the panel motion requires solution for the air pressures as 
an integral part of the calculation. The importance of this 
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coupling has long been recognized as the key factor leading 
to the instability known as panel flutter . In addition, 
some work has suggested that the coupling with a supersonic 
air stream can introduce damping that is large in comparison 
with damping from structural sources. It has also been 
found that coupling of panel motion with a heavy fluid, such 
as water, can cause large reductions of the plate resonant 
frequencies, compared to the in-vacuo natural frequencies.' ' ‘ 
Previous Panel Response Research 

The mathematical difficulties involved in solving the 
equations describing the motion of a flexible plate coupled 
to a compressible, flowing air stream and driven by a statis- 
tically random pressure field have led to the use of approx- 
imations of various kinds. First, the aerodynamic pressure 
field is separated into two parts that do not interact. 

One part is the random pressure field associated with the 
boundary layer turbulence, which is taken to be given from 
experimental measurements. The other part is the additional 
pressures caused by the motion of the panel. The effects of 
interaction of these two pressure fields apparently cannot 
be calculated within the current state of the art. For some 
simple types of structures, such as uniform spherical and 
cylindrical shells and infinite flat plates the panel motion- 
pressure interaction problem is simplified because the fluid 
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pressures don't change or couple the in-vacuo normal modes 
of the structure, or because the panel motion can be 
represented by a simple traveling wave . In some 
studies the panel pressures were not included at all, and 
attention was thus focused on the statistical or plate 
model aspects of the problem. In other studies the ef- 
fects of panel motion pressures were included as an equiva- 
lent damping or radiation resistance factor, and numerical 
values obtained from experiment . ^ Aerodynamic piston 

theory has been used to analyze the coupled response for air 
stream Mach numbers greater than one. This theory sim- 
plifies the coupling of plate motion to aerodynamic pressure 
because it relates the pressure at a point on the plate to 
the motion at only a few points , rather than to the motion at 
all other points of the plate, as is done in the more exact 
aerodynamic theories. ' ' For plates vibrating in water, 
the "virtual mass" feature of the coupling has been included 
but the damping feature not included. Additional research 

on various aspects of this plate response problem is described 
in the papers referenced by Lin, Dowell , Strawderman, 
Maestrello^^®^ and Leibowitz. All the work described in 
the above papers contains the assumption that the gradients ; 
of velocity and temperature associated with the boundary 
layer do not influence the coupled dynamic response Of the 



plate-fluid system. 
Boundary Layer Effects 
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Recent research has shown that the thickness of the 

boundary layer over a panel can have an important effect 

(20 21 ) 

on the flutter properties of the panel. ’ The exper- 

imental results show that the dynamic pressure at flut- 
ter in the presence of a particular boundary layer was twice 
the dynamic pressure value obtained (by extrapolation) for 
zero boundary layer thickness. Calculated results ' also 
show the large influence of boundary layer thickness on 
flutter dynamic pressure. Because of the close connection 
between panel response and panel f lutter , these flutter 
results imply that the characteristics of panel response will 
depend on the thickness of the adjacent boundary layer. Such 
effects might also be expected on intuitive grounds. For 
example, in the case of an extremely thick boundary layer 
the fluid velocities existing at large distances from 
the panel could be expected to have only small effects on the 
panel motion, as compared to a thin boundary layer in which 
the large velocities exist within a fraction of an inch of 
the panel surface. Considerations such as these suggest 
that it would be worthwhile to develop methods for calculating 
the response of panels that are coupled to a passing air 
stream through a boundary layer containing gradients of 
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velocity and temperature. Development of such methods is 

the purpose of the work described in this paper. 

(9 21 ) 

The procedures developed by Dowell ' ' for calculating 
flutter of panels take into account the effects of velocity and 
temperature gradients and are sufficiently general to allow 
calculation of response as well. Dowell's method of solution 
uses a combination of plate modal series expansion, Fourier 
transforms in the direction parallel to the plate and finite 
differences across the boundary layer to reduce the govern- 
ing partial differential equations to ordinary differential 
equations in time which are then solved numerically on a 
digital computer. This procedure is very general and can 
incorporate a great variety of phenomena such as plate non- 
linearity, aerodynamic velocity and temperature gradients, 
and transient input forces. When calculating results for 
particular practical configurations this method is probably 
the best type to use. Such great generality however, tends 
to increase the computation effort required to obtain each 
data point, and consequently increases the difficulty of 
carrying out trend studies to explore the effects of bound- 
ary layer thickness, for example, and increases the difficulty 
of obtaining the statistics of the panel response from the 
statistics of the forcing pressures of the boundary layer. 

An approach alternative to the time domain solution of 
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Dowell is the frequency domain solution. This procedure 
makes use of either a harmonic time dependence or a Fourier 
transform in time instead of the numerical integration in 
time. The frequency domain solution appears to offer the 
advantages of (1) requiring less difficulty in carrying out 
trend studies, (2) allowing more straightforward develop- 
ment of relations between input and output statistics, and 
(3) allowing easier incorporation of simplifications such 
as asymptotic solutions of integrals. The frequency domain 
solution is used in the work described in this paper. 

The Present Research 

The purpose of the work reported herein, then, is to 
develop methods for calculating the response of panels 
coupled to a passing air stream through a boundary layer 
using a frequency domain solution. In order to focus atten- 
tion on the effects of the boundary layer, other aspects of 
the problem are simplified. The structural panel considered 
is uniform, flat, isotropic, elastic, initially un-stressed, 
and responds in the linear, small-deflection, amplitude, 
range. The panel has infinite width and finite length. 

The flow is parallel to the panel length. The applied force 
is harmonic, rather than random. The fluid is considered 
to be compressible and Newtonian; its motion is taken in 
two parts: (1) a large magnitude shear-flow motion parallel 
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to the panel surface, and (2) an acoustic motion of infinit- 
esimal magnitude and arbitrary direction. The large motion 
is taken to be steady, laminar, and unidirectional. Viscos- 
ity effects are included in the large motion, but are not 
included in the small magnitude part of the flow. 

The plan of this paper is: (1) derive appropriate 
partial differential equations and boundary conditions gov- 
erning the plate deflection and the fluid pressure, and 
relate them to previously derived equations; (2) solve the 
governing partial differential equations so that numerical 
values of specific deflections and pressures can be obtained 
using a digital computer; and (3) present results displaying 
the effects of boundary layer thickness on panel deflection 
and aerodynamic pressure. 



Chapter II 

Derivation of the Governing Equations 

In this chapter the equations governing the response 
of a coupled system consisting of a flat plate and a flowing 
fluid are derived. In addition the boundary conditions on 
the fluid, on the plate, and at the plate-fluid interface 
are obtained. A general assumption used is that the response 
is small , so that the equations to be obtained are linear in 
the dependent variables of plate displacement and fluid pres- 
sure. The system under consideration is shown in figure 1. 

A flat rectangular elastic plate occupies a region of the 
x-y plane; the remainder of the x-y plane is occupied by a 
baffle, or rigid plate, that serves to separate the upper 
half space from the lower and provide boundary conditions 
at the edges of the plate and at the lower surface of the 
fluid. A Newtonian flowing fluid fills the upper half space, 
z > 0; the lower half space, z < 0 , is empty. A known force 
q(x,y,t) is applied to the plate surface. This system can 
be thought of as an idealization of, for example, a small 
panel in the surface of an aircraft wing or fuselage. In 
the development of the equations particular attention is 
given to the effects of velocity and temperature gradients 
that form the boundary layer in the fluid near the plate, 
and other aspects of the system are treated in simplifying ways. 
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The Plate Equation 

A sketch of a small element of the plate showing the 

forces acting is presented in figure 2. The equilibriiom of 

(23) 

an element of the plate is governed by the equation : 

4 . . T 

DV w = load intensity !• 

The load intensity acting includes: 

the inertia reaction force, -p — 5 - 

P 3t 

the fluid pressure, -p 

and the given external force, q. 

Inserting these forces on the right hand side of equation 1, 
and rearranging terms leads to the equation governing plate 
motion : 

DV^w + p ^ + p (x,y, 0 ,t) = q (x,y,t) 2 „ 

P 3t^ 

For the present work the pressure p(x,y,0,t) is the pressure 
in the fluid that results from the plate motion w, thus 
to solve equation 2 the relation between p • and w must 
be obtained from the governing fluid equation and must be 
used in the solution. The pressure is a function of posi- 
tion in space and time; the form p (x,y, 0 ,t) indicates that 
the pressure at z = 0 is considered to be applied to the 
plate surface. When there is no fluid present the pressure 
is zero and equation 2 can be solved using known techniques. 



z 



Figure 2.- Plate element and forces. 
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In order to focus attention on the fluid boundary layer 
effects a very simple type of plate has been assumed in ob- 
taining equation 2. The plate is flat, has uniform thickness, 
mass and stiffness distributions, and the plate material is 
isotropic and elastic. There is no initial pre-stress and 
the response is assumed to be in the linear, small amplitude, 
deflection range. This simple plate is felt, however, to 
represent some realistic systems. 

The Aerodynamic Equation 

The equation governing the pressure in the fluid is 
obtained from the basic equations: 

p*D(u*) = p*f^ - Vp* + u* + yV(V*u*)/3 3. 

D(p*) + p* V • u* = 0 4. 


p* = R p* T* 


p* D(T*) 



9T* 

3x 


1 


9u* 




5. 


8u* 




6. 


These equations are, respectively, the Navier-Stokes equa- 
tion governing a Newtonian fluid, the continuity equation, 
the equation of state for an ideal perfect gas , and the 
energy equation for a Newtonian fluid obeying the Fourier 
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law of heat conduction. These six equations (equation 3 is 
a vector equation) are sufficient to determine the six un- 
knowns ; pressure p*, density p*, temperature T*, and 
velocity (vector) u*. All six are necessary to describe 
the behavior of a compressible fluid. 

The unknowns described by equation 3 to 6 are total 
quantities. For example p* (x,y,0,t) is the total pressure 
acting on the plate, whereas the pressure p(x,y,0,t) 
appearing in equation 2 is only that part of the pressure 
caused by the motion of the plate. The unknown "starred" 
quantities are therefore divided into two parts, as follows: 

u* = U + u 

P* = pQ + P 

P* = Pq + P 

T* = T + T 

0 

where the quantities u, p, p, and T are related to the 
motion of the plate, and vanish when there is zero plate 
motion. The quantities U, p^ , p^, and Tq are the values 
that remain when the plate motion disappears; for convenience 
this flow condition is referred to herein as the basic flow. 
Linear equations governing the "plate motion" quantities are 
obtained as follows. Equations 7 are substituted into equa- 
tions 3 to 6 and the terras of each resulting equation are 
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divided into three groups . One group of terms contains only 
basic flow quantities, U, Pq, Pq, Tq, and disappears under 
the assvunption that the basic flow satisfies equations 3 to 
6. Another group of terms contains squares and/or products 
of the "plate motion" quantities, u, p, p, T, and disappears 
under the assumption that u, p, p and T are small enough 
that squares and products are negligible compared to the 
linear terms. The group of terms remaining are the govern- 
ing equations, however two additional simplifications are 
made. First, it is assumed that viscous effects and body 
forces are negligible under conditions of interest, so the 
terms containing viscosity and body forces are deleted. 
(Recall however that viscous effects were retained in the 
equations describing the basic flow.) Second, it is assumed 
that there is insufficient time for any heat flow to take 
place due to plate motion, so terms involving heat flow are 
dropped. The equations resulting from the above steps are; 


Pq + Pq(U * 7u + u • VU) + p d(U) = -Vp 8. 

+ V • (pU + Pq u) = 0 9. 

p = R(PqT + pTq) 10. 

p d(TQ) + + u-VTq) = -PqV-u-pV-U 11. 

where d( ) = + U |iJ- . 


Equations 8 to 11 form a set of six equations that are linear 
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in the six "plate motion" quantities u, p, p, and T; 
therefore these equations together with given distributions 
of the basic flow quantities U, Pq, Pq, Tq and appropriate 
boundary conditions should be sufficient for obtaining a 
solution. 

It should be noted that equation 10 was obtained from 

the equation of state for an ideal gas, and that there is 

(24) 

no comparably simple equation for liquids. In acousti- 

cal studies use is made of a sound speed relation, instead 
of an equation of state, but sound speed relations for 
liquids are not simple. For example, one empirically deter- 
mined polynomial equation for sound speed in water contains 

terms up to fourth power in temperature and up to third 

(25) 

power in pressure. It appears that extension of the 

present research to flowing liquids, such as water, would 
require that special consideration be given to an appropriate 
equation of state. The present work is thus limited to an 
ideal gas as the flowing fluid. 

In order to solve equations 8 to 11 particular dis- 
tributions in space and time must be obtained for the 
basic flow quantities U, p^, p^, and Tq. The 
distributions of interest for the present work are 
those associated with a boundary layer on a flat plate. 
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In addition the basic flow is assumed to be parallel flow; 
that is, the boundary layer thickness is taken to be con- 
stant in the direction of flow rather than increasing along 
the flow direction, as actually occurs. This assumption 
appears reasonable because parallel flow has been found 

/ 2 g ) 

adequate for hydrodynamic stability studies. In 

addition the parallel flow assumption allows significant 
simplification of the equations, while still retaining 
important features of the boundary layer flow field, namely 
the gradients of velocity and temperature in the direction 
normal to the plate. 

Instead of using equations 10 and 11 in further develop 
ments, the alternate equation 

d(p)+u*VpQ+pd(pQ)/pQ = YRTQ|d(p)+u*VpQ+pd(pQ)/pJ 12. 

will be used. The derivation of equation 12 is given in 
Appendix A. The derivation makes use of the continuity, 
state, and energy equations (equations 4, 5, and 6) to 
obtain a pressure-density relation for the total quantities. 
Equations 7 and the procedure described following equations 
7 are then used to obtain equation 12 from this pressure- 
density relation. Equation 12 is used in place of equations 
10 and 11, and clearly displays the presence of gradient 
terms in the pressure-density-sound-speed relation. 
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Equations 8, 9, and 12 form a set of five equations in the 
five unknowns pressure p, density p, and velocity u, and 
will be used in subsequent developments. 

Written out in full, equations 8, 9, and 12 are: 
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9S^ 9lJ P [ 


16. 


9t 


9y ’ 9z “ 9x 

“g^ -gh“ -g>“ -SH -©— 


0 4. p !£o 
p.9t 


8P0 3Po 9Po 

U ^ ^ ^ W 97 

L-®~ Hgh- —g)— . 


= I7 + u |E + V |E + w If t 


--df" ”g>“ 
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9p> 


+ u 


3x 


+ V 


3y 

-Qr 


+ w 


3P( 

3z 


E 


3P( 

-u ^ 
— 


+ E 
Pq 


3pQ 

3p0 


■"<D“ 





17. 


The assumption that the basic flow is steady deletes the under- 
lined terms nirnbered ^ . The assumption that the basic flow 

is parallel implies that V and W are zero and that U depends 
only on z^^^^ therefore the underlined terms numbered(^ are 
deleted. The assumption of parallel flow also implies that 
the mean density , pressure p^ and temperature Tq depend 
only on z, therefore the underlined terms numbered @ are 
deleted. Finally, the usual boundary layer assumption that 
variation of the basic pressure p^ over the boundary layer 
thickness may be neglected leads to the deletion of the 
underlined terms numbered @ . Using the short hand nota- 

tion 


d( ) 


= 1L_) 

at 


+ u 


9( ) 
3x 


the equations resulting from these assumptions are; 


X 1 3p 

d(u) + w ^ 


18. 

d(v) 

- 1 lE 

■ Pq 9y 


19. 

d(w) 

- ^ ^ 
■ P() 32 


20. 

^ 1 j. ,, ^*^0 . „ Fsu . 9v 9w"1 

d(p) + w ^ ^ + _J 

= 0 

21. 



0 


19 


d(p) 


- YRTq 




22 . 


Equations 18-22 can be reduced by simple substitutions, dif- 
ferentiations and application of the operation d ( ) , to the 
single equation for the plate motion pressure p: 


1 

2 

c 


d^(p) 


d(V^p) 


+ 2 


3U 

3Z dxdz 




(|f> = » 


23. 


Equation 23 written out in full is: 


3U 


3U^ ^ „3 £p 

3t9x 9x 


3^ P 

+ 3^ P + 3^ P 

-U 

3 3 3 

3 P + 3.^ p 9 -^p 

9t9x^ 

9t9y^ 9t9z^ 

9x^ 9x9y^ 9x9z^ 


+ 2 H 1 !^0 

9z 9x9z 


.2 .2 
i. P..... + u ^ P 
9t9z ^ 9x9z 


= 0 


24, 


Equation 24 is the equation governing the plate motion 

pressure p, and is the equation to be used throughout the 

remainder of this paper. The basic flow quantities required 

for solution of equation 24 are the sound speed c, velocity 

U, and density Pq, all of which are functions of the space 

coordinate z as is appropriate for a compressible mean flow. 

f 21) 

Equation 23 has been obtained previously by Dowell, 

and for an incompressible basic flow equation 24 reduces to 

. f 27 ) 

the equation obtained by Graham and Graham. ' ’ 


For classical 
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inviscid aerodynamic basic flow the gradients of U, Pq and 
Tq vanish and equation 18 to 22 can be shown to reduce to 
the classical convected wave equation: 


2U 

3t2 


3x3t 



2 


c 





25. 


Boundary and Initial Conditions 

The solutions sought are for "steady" acoustic behavior 
in a steady basic flow. The steady state is considered to 
have been established by starting from complete rest through- 
out the system and continuing until transients have all died 
out. The boundary of the fluid is taken to be so far away 
from the plate that disturbances originating at the plate 
never reach the boundary : consequently one of the space 

boundary conditions is that for fluid locations far from 
the plate the acoustic disturbances must have only outgoing 
wave parts, or must be decreasing with increasing distance. 

The boundary condition at the surface where the fluid contacts 
the plate is basically that the fluid particles don't pene- 
trate the plate. For the solution of the basic flow problem 
where viscosity has been retained this implies the "zero 
slip" boundary condition ordinarily used to derive boundary 
layer flow solutions. For the acoustic motion the boundary 
condition used allows particle velocities that are 
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tangent to the deflected plate surface. This is done to 
avoid the complications associated with boundary layer so- 
lutions for the acoustic flow. The fluid velocity at any 
point is given by 

u* = U + u, 

= (U + u) i + V 2 + wk. 

Resolving velocities at the plate surfaces in figure 3 and 
requiring the velocity of the fluid normal to the plate to 


Z 



Figure 3.- Sketch of the velocity components at the plate 
surface . 


be equal to the component of plate motion normal to plate, 
leads to 

9w 

w cos 6 ~[U(z) + u] sin 6 = cos 0 26. 


where w = plate displacement. But, linearizing implies 
P 

0, w , u, w << 1 and 

P 9w 

cos 0=1, sin 0 = tan 0 = 


27. 


This condition reduces equation 26 to 

9w 9w 

w| -(0(z) + u] I ^ ^ 

z=0 z=0 


28 . 
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3w 

But u - 5 —^ is a product of small quantities, and U(z) | = 

oX Z“U 

0. 

Finally the boundary condition is 

(x,y,t) 

w(x,y,z,t) I = 29. 

z =0 


To introduce pressure in place of w, introduce the momentum 
equation : 


3w .. ^ 1 

3t 3x Pq 3z 


30. 


On the plate surface U = 0, so equation 30 becomes 

3w _ _ 1 0 l£ 

3t “ Pq 3z 


Combining equations 29 and 31 yields 


3w 

3t 


z =0 



i I 

1^0 2=0 


32. 


This equation supplies the connection between the homogeneous 
equation governing fluid behavior (pressure) and the non- 
homogeneous plate equation. 

For the case of constant basic flow, equation 28 becomes 

3w dw 

«(x,y, 2 ,t) 1^^^ = ^ + U ^ + U |j) Wp. 33. 

Combining equations 30 and 33 leads to 
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I# - (It- + u|^) W - + Ug^) (g^ + Ug^) W^. 


Pq 3z 


'3t 


So the boundary condition for the constant basic velocity case 
is 



PqCO) 


3t 




3t 




w 


34o 



Chapter III 

Solutions for Plate Response and Acoustic Pressure 


In this chapter solutions are obtained for the system 
of equations and boundary conditions obtained in Chapter II. 
The solutions obtained are also shown to reduce to known 
special cases. Attention is focused on the effects of the 
boundary layer; other aspects of the solution are simplified. 
For continuity of presentation the plate equation solution 
and the acoustic equation solution are each given as a 
separate unit, however, it will be observed that the solu- 
tions of the two equations are interdependent. Thus the 
effect of the fluid pressure due to plate motion is inclu- 
ded in the solution for plate motion. Therefore, these 
solutions could be used to study panel flutter or stability 
of a boundary layer as well as panel response. 

The Governing Non-Dimensional Equations 

The equations to be solved are given in Chapter II, 
equation 2 governs plate motion, equation 24 governs the 
fluid pressure, and equation 32 governs the boundary condi- 
tion that couples the plate motion to the acoustic pressure. 
Additional boundary conditions require that waves be out- 
going at large distance above the plate, and that plate 
displacement be zero in the region of the x-y plane occupied 
by the baffle (see fig. 1) . For simplicity only harmonic 


24 
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forces, and therefore harmonic responses, will be considered. 
Only the steady state part of the solution, not the transient 
part, will be determined. In addition, the plate will be 
taken to be infinitely wide and the response will be taken to 
be independent of the y coordinate. These assumptions are 
expressed mathematically by putting the force in the form; 

q(x,y,t) = q(x) exp(iajt) , 35. 

and by putting the responses in the form 

w (x,y,t) = w(x) exp(io)t) 
p(x,y,z,t) = p(x,z) exp(io3t) 

Greater efficiency in the solution will be 
the following non-dimensional quantities; 

A 

X = XL 

A 

z = zL 

A 

w = wL 

p = Pq (0) w^L^p 

o2 _ 2 4 

^ - PpCO L^/D 

y = Pq(0)L/Pp 

q = qL^/D 

Note that in these equations D is the plate stiffness pa- 
rameter (not the differential operator) and y is the mass 


36. 


obtained by using 


37, 
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ratio (not viscosity) « Substitution of equations 35, 36, 
and 37 into the dimensional governing equations leads to 
the non-dimensional governing equations: 

A'' 

rv ** ^ ^ ^ y\ ys 

+ yfi'* p(x,0) = q 38 

3x^ 


3z 2=0 


= w (on the plate) 


39 


and 


/s'' 2"' 3 

■ip -(3U/Lw)-^ + 3i (U/Lw)^ + (U/Lm)^ 


3x 


3x 


3x^ J 


i (c/Lw)^ 

'iS + iS' 

- (U/Lw) (c/Lw) ^ 

r 3"' 3'' 1 


_3x"^ 3z^_ 


_3x"^ 3x3z"^_ 


"" ni ^ ^ ^ p 

+2 (C/Lw)^ 3(U/Lco)/3z + (c/Lw) ^ 

L J 3x3z ^0 3z 


i 

'0 3z L 3z 


(U/Lco) 


2^ 

3 P 

A * A 

3x3z 


= 0 


40 


In the following sections of this Chapter solutions are de- 
veloped, first for the plate equation 38, then for the 
fluid acoustic equation 40 using boundary condition 39. 
Solution of the Plate Equation 

Reduction to Matrix Equation .- The plate differential 
equation 38 can be reduced to a matrix equation by means of 
an eigenfunction expansion as follows. First orthogonal 
eigenfunctions are obtained for the differential operator 
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by solution of 


4'' 

3 /s 

ni , , 

. 11 ... ^ 1 . - — 4 }ifj 

» 4 mm 

3 X 


41, 


Simple support boundary conditions are used herein, 
The orthogonality condition requires that 


/V A. /\ /N 


w (x) w (x) = 1 when m = n 
n m 

= 0 when m ^ n 


42 


/N /V 


Then the plate displacement w (x) is expressed in the series 


A A 


of eigenfunctions w^ (x) as 


w (x) = 


N 


A A 


a^ w (x) 
n n 


n=l 


where the unknown coefficients become the dependent 

A A 

variables to be found instead of w(x). Substitution of 
equation 43 into equation 38, followed by the standard use 
of the orthogonality condition 42 leads to the equation; 


A A 


^m^^m " J p(x,0) w^(x) dx = 


44, 


where the generalized force is given by 


A A A A 


Qm = Iq 


45o 


Further reduction of equation 44 can be made using knowledge 

A A 

of the form of the solution of pressure p(x,0)» As shown 
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in the next section of this Chapter, solution of the acoustic 

/N /S 

equation for p(x,0) can be carried out using the Fourier 
transform pair 


O ^ ^ f 

p(k,z) = J 


/V A 


p(x,z) exp [-ikx] dx 


a 


46. 


00 

^ ^ ^ *1 t Q ^ ^ 

and p(x,z) I p(k,z) exp [ikx] 


dk 


and the series expression for transformed pressure p 

N 

O A • ' ^ Q A A 

P(k,z) =2_^ a^ Pj^(k,z) 


47, 


n=l 


The a^ in equation 47 are the same as the a^ in equation 43. 
Substitution of expressions 46b and 47 into equation 44, 

interchange of the order of integration and summation, and 
use of the definitions 


A A 


Zm(k) =J exp (ikx) dx 


and 


B 


mn 


00 


Pj, (k,0) dk 


48. 


49 


reduces equation 44 to the form 


N 


a (X “ ^ B a = 

m m mn n 


m 


n=l 


50, 


Equation 50 is the basic equation to be solved for the un- 
known coefficients a^. The a determine the plate 

in m ^ 
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displacement and the fluid pressures through the relations 

N 

w(x,y,t) = L a^ (x) 51. 


/V /\ 


p(x,y,z,t) = 


n=l 


(Pq(0)AV2tt) e^“^ 


N «> 


n=l 


ikx ,r 
z) e dk 


52. 


Equation 51 has been obtained by combination of equations 

✓V /V 

36, 37, and 43; the must be found from the eigen- 

value problem equation 41 together with boundary conditions 
at the edges of the plate. Equation 52 has been obtained 
by combination of equations 36, 37, 46, and 47; the functions 

O /N A 

Pj^(k,z) must of course be obtained by solution of the acous- 
tic equation. 

When the parameters entering the dimensionless equations 
38 and 40 are given, the only unknowns in equation 50 are 
the a^^. Then equation 50 can be used to form a matrix equa- 
tion that can be solved for the vector {a^ a 2 ...ajj) using a 
procedure such as Gauss elimination. Equation 50 contains 
the coupling between the plate motion and the acoustic (aero- 
dynamic) pressure in the term , therefore the stability 

properties (flutter, hydrodynamic instability) of the plate- 
fluid system can be obtained by appropriate solution proce- 
dures. In the present study, however, the system is assumed 
to be stable and the solution sought is the magnitude 
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of the plate displacement and the acoustic pressure when 
a given force is applied. In addition, the present 

study is concerned with the effects of boundary layer 
profiles on response, rather than with the complexities 
of solution of equation 50, therefore approximate response 
solutions of 50 are developed in the following sections. 

Solution for small mass ratio y .- Equation 50 holds 
for each value of m from 1 to N. When all N equations are 
written out they can be arranged in the following matrix 
form: 

["a-]{a} + y [B] {a} = {Q} , 53. 



To solve equation 53 it is assvuned that y << 1, so the 
solution can be expressed as 


{a} = {a>Q + y{a}j_ + y^{a >2 + 


55. 
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Substituting 55 into 53 and equating coefficients of each 
power of y leads to the solution 

{a} = [a]"^ ([I] - y[B][a]"^) {Q} + O(y^) 56„ 


When higher order terms are dropped^ equation 56 yields 
the following expression for each a^. 



N 



B 


mn 






a^) 1 


n=l 


V< 


m 




57, 


Using the B^^^ found from solution of the acoustic equation^ 
the a^ from equation 57 are used in equations 51 and 52 to 
obtain the plate displacement and the acoustic pressure. 
Solution of the Acoustic (Aerodynamic) Equation 

The acoustic pressure is governed by equation 40, sub- 
ject to the boundary condition equation 39 and a radiation 
condition. The procedure for solution of equation 40 uses 
a Fourier transform in the x variable to reduce the partial 
differential equation to an ordinary differential equation 
whose solutions are expressed in transfer matrix form. 

Reduction to an Ordinary Differential Equation .- The 
transform pair for pressure, equation 46, reduces the part- 
ial differential equation 40 as follows. Equation 40 is 

yv A 

multiplied by exp(-ikx) d5^, integrated from - to «> , and 
the integrals reduced by integration by parts. The resulting 
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o 

equation governing the transformed pressure p is 


- 


2Mc 

— X — 


1+cMU 3z 


3U 1 ^^0 

P 


'0 3z 


3z 


\ 2 

1+McU 


/V 

c c 



p = 0 


58. 


Similar transformation of the boundary condition equation 39 , 

/V 

and use of the series expression for w(x) yields 


O 


^ I 

3z 


z=0 



n=l 


(k) a 

n 


59. 


Equation 58 is linear and homogeneous, therefore the solution 
will depend linearly on the boundary conditions, equation 
59. As a result a solution can be obtained for each term of 
the boundary condition series and the complete solution is 
then the sum of these separate solutions. For convenience 
the a^ are included in the series expression as follows : 


p(k,z) 



n=l 


O 



(k,z) 


(47.) 


It is now clear that each p^ must satisfy equation 58 and 
the boundary condition 



3z 


z=0 


Z* 

n 


(k) 


60. 


Additional radiation boundary conditions must also be inclu- 

O 

ded. We now observe that when p^ is obtained from equations 
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58 and 60, the result is to be inserted into equation 49 to 

obtain the B which are then used to obtain the a coef- 
mn 

ficients from equation 57. Finally the a^^ determine the 
plate displacement and pressure by equations 51 and 52. 

General Solution by Transfer Matrix .- Equation 58 is 
an ordinary differential equation with variable coefficients. 
The coefficients are variable as a result of variations, 
through the boundary layer, of velocity U = U(2)/U^, density 
Pq(z), and sound speed c = c(z)/c^. In general the solu- 
tions will not be obtainable as known functions, so numerical 
procedures must be used. In addition the boundary conditions 
are in two parts, that is, to obtain p at z=0 one boundary 
condition on 3p/9z is given at z=0 and the radiation condi- 
tion is applied at some large value of z. A convenient means 
of expressing the solution for these conditions is transfer 
matrices. Before carrying out the solution the radia- 

tion condition will be developed. 

For the profiles of interest in this study the coef- 
ficients of equation 58 are constant for very large values 
of z because of the boundary layer nature of the basic flow 
^ velocity, density and sound speed. In order to apply the 
radiation condition it is assumed that U, c and are 
exactly constant for all values of z > 6 . Then 6 is the 

dimensionless boundary layer thickness. It is now convenient 

o ^ o 

to change the independent variable to z where: z = o z . 
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The governing equation and boundary conditions are now: 



/N 

2Mc 

3U ^ 3Pq 

o 

3p 

+ 

'/i+mJu\ 2 1 

1+McU 

3z ^0 32 

3z 

1 

1 

1 u 

C O 


ok p 


n 


0 , 


61. 


and 


o 



3z 


o 

z=0 


6 Z*(k) 


62. 


o 

Now, outside the boundary layer (z > 1) equation 61 reduces 
to 


2 ° 

9 P 

~~g~ 9 ' + 

3z^ 


1 2 
( 4 - + M) ^ - 1 


0 ''0 o 

<S k'^ = 0 


63. 


Equation 63 has constant coefficients and solutions 


exp (iaz) + exp (-iaz) , 


If 


= 


(1 + M )^ -1 


2 "'2 

6 ^ > 0 ; 


and 


n 


o o 

Cj^ exp (bz) + exp (-bz) , 


If 


2 2 
a^ = -b*^ < 0. 


64. 


65. 


66 , 


67. 


Substitution of solution equation 64 into the pressure ex- 
pression equation 52 shows that only outgoing waves will be 
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present if = 0. Outgoing waves have the form exp (icot - 

O 

iaz). Substitution of equation 66 into the pressure equation 
52 shows that only decreasing pressures will be present if 
= 0. Now the procedure for solution of equation 61 sub- 
ject to the "radiation" condition and boundary condition 
equation 62 can be summarized as follows. Equation 61 will 

o o 

be solved for 0 < z < 1, the resulting p^ will be forced to 

O 

satisfy boundary condition equation 62 at z=0 , and be matched 
o 

at z = 1 to the solution 


p^ (k,z) = exp (-iaz) , 

if a^ > 0 ; 

or to the solution 


if 


O o o 

p^ (k,z) = exp (-bz) , 


2 

a 


< 


0 . 


68 . 


69. 


This procedure can be carried out efficiently using the trans- 
fer matrix technique. The transfer matrix idea is to relate 

O o 

the pressure p^ and. the derivative of pressure p' at any 
o 

two values of z by the relation 



The values of the T^j depend of course on the two values of 
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z. To apply the boundary conditions at z^= 1 and Z 2 = 0 
we use equations 68, 69 and 62 in 70 to obtain 


B e 
n 


-xa 


-xa 


- 


O 

T T 

11 12 



T„, T^.^ 


6 Z* 

21 22 


n 


71, 


and 


when a > 0 ; 




-b 


1 

-b 


m m 

11 12 

T T 

21 22 


o ^ 




6 Z* 
n 


72, 


when a <0. 

Equations 71 and 72 each form a set of two linear algebraic 

o o 

equations in the two unknowns p and B , and p and D , 

n n n n 

o 

respectively, and can be solved for p^ to obtain: 


Pn (k,0) = -6Z* (ia ’’ll + ’'21> ' 


73, 


when a > 0 ; 


and 


(k,0) = 


-«ZJ (b + T22>/(b 


74 


2 

when a <0. 

Derivation of the for given coefficients of equation 61 
can be carried out using standard numerical techniques; 
examples are presented in later sections of this paper. The 

O 

expressions 73 and 74 for Pj^are now the final results required 
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for the computation of plate displacement and pressure 
using equations 51 and 52. 

O '' 

The solution procedure is as follows. Values for p^(k,0) 

from equations 73 and 74 are inserted into equation 49 to- 

gather with the mode transform Z (k) to obtain values of B, . 

^ m mn 

These values are used in equation 57, along with the general- 
ized forces and values of the parameters, to obtain the 
series coefficients a . The a are then used in equations 
51 and 52 to obtain the plate displacement and acoustic pres- 
sure. In order to obtain pressure at a non- zero value of z 


O /V 

from equation 52, values of p^(k,z) are obtained from the 

O 

transfer matrix equation 70 using the particular values z^^ = 
^ o ^ ^ o o ^ 

z/6 to obtain p (k,z) , and z, = 0 where p (k,0)is known. 

Solution for T^^^ by Numerical Integration .- In general 
the coefficients of equation 61 will be sufficiently comp- 
licated that analytical or series solutions will be imprac- 
tical, therefore numerical integration must be considered. 
Two regions of numerical solution must be considered; one 
region near the singular point of equation 61, and the other 

o 

region covering the rest of the interval 0 < z < 1. 

The singularity of equation 61 is associated with the 
zero of the denominator of the coefficient of the first de- 

/V 

rivative term. The location of this singularity on the k 


axis is of interest because of the integrations on k that 

must be carried out to obtain B and pressure. The 

mn ^ 


38 


o o 

coefficient of 9p^/9z is singular where 


1 + MXk U (z) = 0 


75 


In this expression c = kX has been used. Equation 75 shows 

/S A. /N 

that when k > 0 there is no singularity; when k < 0 but |k| 

_ O 

< 1/(XM) there is no singularity because the velocity U(z) 

_ o ^ 

has values in the interval 0 < U(z) <1; and when k < 0 and 
/\ o 

|k| > 1/ (XM) the singularity is located at values of z 

obtained from the equation 


U(z) = l/(MX/|k| 


76 


A sketch of these results is shown in figure 4. 



Figure 4.- Location of singularity in equation 61. 
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The particular shape of the curve locating the singularity 

_ o 

will depend on the particular variation of U(z). Solution 
of equation 61 in the region near the singularity could be 
carried out using a Taylor series expansion of the coeffici- 
ents and a power series solution of the equation. 

Solution of equation 61 in the region away from the 
singularity can be carried out using the procedure of ref- 
erence 28. To begin, change notation so that equation 61 
appears as 

2 

— - a(z) — ^ + b (z) p = 0. 77. 

3z^ 3z 

o o 

The definition of a(z) , b(z) and p can be seen by compari- 
son of equation 77 witn equation 61. wext define 

3d ® 

= p , Z 2 = -fr , and z = s 78. 

3 z 


so that now equations 77 and the second of these relations 
can be written 


3 

3s 




79. 


Equation 79 is the same form as equation 5-45 of reference 
28, therefore 79 has the solution presented in the reference. 
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Writing 

{Z} 



80. 


the solution obtained from reference 28 is 

^n+1 = ^n I <^0 + 2k2 + k3> 

Where k„=hA(s„)z„ = hA (s„ ^ |)(z^ + 

’=2 = hA (s„ + |)(Z„ + 

kj * hA (s„ + h)(z_^ + k,) 

r A(Sjj + h) , and A(Sj^ + h/2) are, respectively, the 
values of the matrix A(s) at the sections n and n + 1 and 
at the section midway between n and n + 1. 

Combining the above unnumbered equations gives the 
relation between z . , and z 

■riAl ^ • 


'n+1 


= < I + 


|[a(s„) + 4a(s„ + I) + A(S„ + h)] 

1“ [* (=n * I) *<V + '^l^n *>>*(®n 

(=n + I)] I2 ( ^n + I) + 





41 


+ A(s^ + h) 



A(s^ + h) X 


A 


2 




Z 


n 


81. 


where the contents of the braces represents the transfer 
matrix linking the state vectors and Equation 81 

relates Z at two closely spaced points. To relate Z at 

O 

two arbitrarily spaced points in z, write 81 as 




and then combine by substitution for several values of n 
to obtain 




[t2] [tj_] [tQ]{Z}Q, 


82, 


where N is the number of the station corresponding to the 

O 

desired value of z and 0 denotes the station where the 
vector {Z} is known. Reverting to the notation 78 and 
combining the product of the [t^ ] matrices into a single 
matrix, equation 82 appears as 


p 

% 

1 

^11 ^12 
T T 


P 

dZ 

lo 0 

21 ^22 


8 z )o Q 


Z = Zj^ 



z=z 


83. 


The Tj^j appearing in equation 83 are the required transfer 
matrix elements for use in equation 70. 
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Solution for T. . by Series for Mach No. < 1.- Tempera- 


ture profiles presented in reference 29 indicate that the 
rise of temperature in a compressible boundary layer is not 
large when Mach number is less than 1. For these Mach num- 
bers a reasonable approximation is that sound speed ratio, 
c = c(z)/c^ , is unity and that density is constant. In 

addition, the velocity profiles presented in reference 29 
suggest that the simplest representation having the charac- 
ter of a boundary layer profile is: 

_ o o 

velocity ratio U = U(z)/U^ = z 84. 

While this representation is not a closely accurate descrip- 
tion, it does have two advantageous uses. First, equation 
84 does represent the general character of the boundary 
layer profile, therefore results obtained using it are ex- 
pected to be at least a first approximation to the correct 

results. This linear description has been used previously 

(27) 

in acoustic studies. Second, use of the description 

84 leads to a reasonably simple equation governing pressure, 
such that solutions can be obtained in series form. These 
series results can be compared with results from the numer- 
ical integration procedure to help gain confidence in the 
results of both procedures by providing two independent 
solutions. Having this confidence the numerical integration 
procedure can then be used to obtain more precise results 
using more realistic representation of the boundary layer 
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profile. No essential increase of difficulty appears to be 
associated with the inclusion of temperature variations in 
the numerical integration technique. 

Using the representation c = 1, Pq “ constant and U = 
equation 61 becomes 


3 p 
°2 


n 


2Mc ^Pn . 
— TTtr —c + 


1+Mcz 3z 


^o\ 

1+Mcz 


Introducing the notation 


- 1 


"2:2 ® . 

c 6 Pn = 0 


85, 


X = 1 + Mcz 


^2 

a = 6 /M c 


o 

P = P 


n 


''2 2 
8 = 6 /M^ 


86 , 


reduces 85 to 


p" - 2xp' - 


(Sx^ - ax'*) p = 0 


87, 


Equation 87 has a regular singular point at x = 0 , and can 
be solved in series form in the neighborhood of the singular- 
ity by the method described in reference 30, Theorem 4.4, 
p. 188 ( Frobenius method) . The solution obtained is: 


p(x) = A pj^ (x) + B P 2 (x) 


88 . 


where 


p^ (x) = I x| j 1 + a. 


n=4 


a^ xV[n (n+3) ] | 


89, 


n even 


N O 
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with aQ = 1, 32 = .le, = 3 a ^_2 - a 


and 


P2(x) = 1 + a2 ^ ^ xV[n(n - 3)] 

n=4 

n even 

with aQ = 1, 

a2 = -.53, 

a = 3a „-aa * 
n n-2 n-4 


90 


From the theorem of reference 30, the series appearing in 
equations 89 and 90 are convergent. Equations 89 and 90 
show that Pj^ (x) and ?2 (^) are defined for both positive and 
negative values of x, and are finite and continuous at x = 
0. Because Mach number M appears in the denominator of a 
and 3, the solution 89 and 90 may converge slowly for small 
M. The reduction of this solution to the limits M = 0 or 
6 = 0 , M ^ 0 to demonstrate reduction to known solution 
does not appear to be straightforward. For these reasons 
a second series solution is developed. 

In the neighborhood of an ordinary point of equation 
85 a series solution can be obtained by the method of ref- 
erence 30, theorem 4.1. The solution obtained by that 
method is ; 

p^ (z) = c f^ (z) + D f 2 (z) 91 


where 
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f ( 2 ) 

= 1 + a 

with 

a 


a 


a 


°2 

z + a. 


°3 ^ 
z + a, 


°4 ^ 
z + 


00 

E 

n=5 


a z' 
n 




= 1 = 0 
= -i 6^[1 - 

^ o r A < 


and 


= I 


;^2r, ;^2,2 ^2 ^2 

01 . 1 -c] -M c 


^ , 0 . ° °2 °3 °4 V '* On 

£ 2 ( 2 ) = z + a 2 2 + 2 + a^ z + ^ ^ a^ z 

n=5 


with 


ao = 0, 


a2 = Me 


^ 1 
3 


1 


„2 "2 1 P n "2, 

M c “J'S [1-c] 


^4 = Mc6^ [2 - c^] 
and, for both f^^(z) and £ 2 ( 2 ) 


^n = 


92, 


93 


i A ^ 

-(n-1) (n-2)Mc + 2(n-l)Mc a^_^ -[1-c^] 6^ a ^^_2 

-[3-c^] Mc6 ^ a ^_3 -3M^J^ 6^ a^_^ - P 6^ 
/[n(n-l) ]. 


n“5 


The function Pj^( 2 ) defined by equations 91-93 can be used 
to obtain the transfer matrices for use in equation 70. The 
particular transfer matrix relating pressure on the plate 

O 

surface z = 0 to pressure at the ddge of the boundary layer 
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o o 

z = 1 is obtained as follows. The relation for and 

its derivative can be written as 


O 




f^iz) 

f[(h 


f2(z) 

f^(z) 



94 . 


O 

Evaluation of the f matrix for z = 0 using equations 92 
and 93 shows that 



95. 


Combining equations 94 and 95 and evaluating at z = 1 shows 
that the transfer matrix elements required for solution of 
equations 73 and 74 are given by 

96. 


T 

■^11 

T 

12 


f^d) 

1 

rH 

CM 

4.1 

T 

21 

T 

22 


f[{l) 

f^d) 


Solution for Zero Boundary Layer Thickness 

The solution procedures developed in this Chapter can 
be checked by comparing a pressure distribution for zero 
boundary layer thickness with the pressure distribution pre 
sented in reference 31 based on inviscid aerodynamic flow. 
The pressures given in reference 31 are based on an elastic 
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plate of infinite extent undergoing displacements in the 
form of a traveling harmonic wave. The pressures for this 
situation are obtained using the method of this thesis 
starting with equations 39 and 40. Equations 39 and 40, and 
the "radiation" boundary condition, are sufficient to deter- 

/s 

mine the pressures for a given plate displacement w . A 
traveling wave in the plate can be represented as 

w (x,y, t) = a^ exp i (cot + kx) , 97. 

and the pressure solution corresponding to this plate motion 
can be written as 

p (x,y,z,t) = c(z) exp i (wt + kx) 98. 

Converting to non-dimensional quantities and substituting 
equations 97 and 98 into equations 39 and 40 leads to the 
results that a(z) must satisfy the ordinary differential 
equation 58, (and therefore 85) , subject to the boundary 
conditions of "radiation" at z -*■ ® and 


9a 


■"~yr 


9z 


z=0 


99. 


This boundary value problem is the same as that satisfied 

O '' 

when 2*(k) is replaced by "^o"' seen 

from equations 58-60. (This result suggests the following 

O Q O 0^0 

interpretation of the p^(k,z): p^(k,z) represents the 
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pressure in the fluid due to an infinitely long harmonic 
wave on the plate surface traveling with wave speed w/k 

/V 

and wavelength 1/k. The function Z*(k) determines the 
amplitude of each wave component so that the svim of all the 
traveling waves on the plate will form the model function 

/S A O 

w^(x).) Therefore the solution derived for p^(k,0)can be 
used for a, with Z* replaced by "ag"* From equation 74, the 

A 

solution for a is 


a(k,0) = -6ag (b "^ll '^21^ 

2 

Note that this solution, derived for a <0, reduces to 

2 

the solution (equation 73) derived previously for a >0. 

The T. . are related to the series solutions of the governing 


O O 


equation 58, 


and 

f2(z) 

by 




^11 

1 

(N 

iH 



f^d) 

f2(l) 


'^21 

T 

22 



f[{l) 

f^(l)_ 


The function and their derivatives can be evaluated 

A 

for zero boundary layer thickness 6=0 from equations 92 and 
93 to obtain 


T T 

•^11 12 


1 2 ^^2 

1 1 + Me + i M c 



'' 2^2 

T T 

21 22 


0 1 +2Mc + M c 

* 


- 


101. 
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Substitution of the from equation 101 into equation 100 

yields 

a{ic,0) = -X Sq (1+Mc)/'/ c^- (1+Mc)^ 102. 

Substitution of 102 into the dimensional pressure expression, 
equation 98, yields 

p(x,0,t) = ~Pq w^L^X aQ(l+Mc)^/ / c^-(l+Mc)^ exp i (wt+kx) 

103. 

Comparison of this result with the pressure expression given 
in reference 31 must take account of the use by reference 31 
of a coordinate system moving at the free stream velocity U; 
and the use of a wave moving in the positive x direction, 
whereas the wave described by equation 97 is moving in the 
negative x direction. The coordinate relative to the moving 
axis x' is related to the fixed axis x by 

x' = X - Ut 

and the speed of the wave relative to x' is 

C - U 104. 

/s 

Using the definition of c leads to 

<-Mc. 

so the pressure expression equation 103 reduces to 

_1 

p(x,0,t) = -Poc'^ Sq k [l-^|— ^ ^ exp i(-k) jc't-x'j 


106. 
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When the names of the variables in equation 106 are changed 
to agree with the notation used in reference 31, the pres- 
sure expression, equation 106, is the same as the expression 
for surface pressure obtained by combining equations (2.4a) 
and (3.5) of reference 31. This result lends confidence in 
the correctness of the analysis methods used in this paper. 
Computations 

The computation procedures used to obtain numerical 
values of plate displacement and acoustic pressure are de- 
scribed in this section. The primary quantities required 
to obtain displacement from equation 51 and pressure from 
equation 52 are the coefficients a^ and the pressure 

0 ^ A 

quantities p^ (k , z ) . 

The coefficients a are determined from equation 57, 

m 

the approximation for small p ; 


m 


1 - y 


N 

■E 


B /[Q (X -n^) 3 

mn n '• m n 


n=l 


Q /(X -n' 
m' m 


(57) 


To determine a the quantities B and Q are needed: 
m mn m 


00 



II 

Pn(K.O) Z^(k)dk, 

(49) 


— 00 



where 

= f 

/s 

'' ikx 

Wm(x) e^ *dx; 

(48) 


0 


‘m 


I 


q(x) Wj^(x) dx. 


and 


(45) 
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The eigenfunctions w^(x) and eigenvalues are obtained 

by solution of 


9 (x) 


= X w (x) 
mm 


subject to boundary conditions at x = 0 and x = 1. For the 
conventional simply supported boundary conditions where 
plate deflection and moment are zero at both ends of the 
plate the solutions of 41 are 


w (x) = /2 sin m 7T X, X = (mu) 
m m 


From these eigenfunctions equation 48 yields 


2„(i^) = mu ^ -i 


2 2 r2 

m IT -k 


2 2 T2 
m TT -k 


The real and imaginary parts of Zj^(k) are illustrated in 

figure 5. Figure 5 shows that the Z functions decrease 

m 

/N 

in magnitude quickly as k increases, indicating that inte- 
grals such as equation 49 need not necessarily cover a 

/s 

large range of k. 

The force distribution used in these computations is 


q(x) = /2 sin ttx 


so the generalized force Q^, equation 45, is 


= 1 for m = 1 
m 


Qn, = 0 
m 


m 1 
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Figure 5.- Transform of the mode function. 
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Using this result, and the eigenvalues from equation 107, 
in equation 57 leads to the expressions for the first few 
coefficients a 

m 


a^ = 1 - U n^/(TT^ - 

/(TT^ - 

= -y ^ / 

4 2 4 4 2 

(tt^ - (2^tt^ - 

^2 , 

a^ = -y B^j^ / 

/ ^ ri2\ ^ n2\ 

(tt - ) (3 tt - fi ) 


a 

b 


c 


111 . 


2 4 

For n very nearly equal to it the first coefficient a^^ is 
large compared to the other a^^, therefore for this study 
it is assumed that a^^ = 0 for m ^ 1 and equations 111a is 
used for a^. Using this result the plate displacement and 
pressure are written: 


/o • iwt 

w(x,t) = a^/2 sin ttx e 


112 . 


• ? o - o • - 

p(x,z,t) =^Pq (0)a)^L^/2Trj a^^ J dK 


The next quantities required are the p^(k,z) so the integrals 
in equations 112 and 49 (the integral required for cal- 
culation of a^) can be evaluated. The expressions for the 

O ^ O 

p^(k,0) are (only p^^ is needed); 

p^(ic,0) = -6 Z* (ia Tj ^2 + *^22^/^^^ "^ll '^21^ 

2 

when a > 0 ; 
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p^{k,0) = -6 Z* (b '^11 '^21^ 

when a^ < 0 . 

The transfer elements are given by 


11 

T 

12 

II 

l-h 

I-* 

f2(l) 

21 

T 

22 

f^d) 

f-(l) 


and the functions f^ and ±2 used here are the series expres- 
sions derived for small Mach number, equations 92 and 93. 

O A /\ 

To evaluate p^(k,z) the general form of the transfer matrix, 
equation 70, is used, to obtain 

p^(k,z) = Tj^^ p^(k,0) + T ^2 

For these calculations pressures were calculated only on the 

o o 

plate z=0 and at the edge of the boundary layer z=l. On the 


plate surface 


’’ll = 1 ’'12 = “ 


and at the edge of the boundary layer 


'^11 "^ ll ' "^12 "^12 


The real and imaginary parts of p^(k,0) are defined by 


p^(k,0) = -X 


6p 6p 

'1 [~5r ’ -JT_ 


The quantities 


are shown in figures 6 and 7 
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plotted as a function of c for Mach number of 0 . 5 and sev- 
eral values of 6/X. Figure 6 shows that for moderate values 

<SPc 

of 6/X the curves for -j — are smooth, so no numerical dif- 
ficulties were anticipated in numerical integration of this 
function (and none were encountered) . Figure 6 also shows 
that, as the value of 6/X becomes smaller, the curves ap- 
proach the curve for 6/X = 0 . This result provides numerical 
support for the result obtained previously by analytical me- 
thods . 

6p 

Figure 7 shows that the real part is a smooth func- 

✓v 

tion except in the region c - -.65. Further examination 
6p o ^ 

shows that -y-/ and therefore pj^(k,0), has a singularity in 

A A 

this region of c; the exact value of c depends on the Mach 

o ^ 

number and on 6/X. Examination of the equation for Pj^(k,0), 
equation 74, shows that the singularity is associated with 
a zero of the denominator, i.e. 

b = 0, 117. 

and is of the form | c - Cq ( . Numerical values of the inte- 

o ^ 

grands of and 1^^ (an integral involving p^(k,0) to be 
discussed in a later section of this paper) are shown in 
figure 8. Figure 8 indicates also that the singularity has 
the form |c - Cq| . The appearance of this singularity for 

A 

a real value of c indicates that the function being trans- 

A A 

formed - p(x,0) - has some sort of non-regular behavior, such 
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as not becoming zero for large x, for the value of c at 
the singularity. Consequently the transforms, equation 46, 
must be regarded in a generalized sense and their existence 
established by an extension of the Fourier integral theorem 
as discussed in reference 33. The result of this general- 

yv A 

ization is that the variable c or k must be regarded as a 

o ^ 

complex variable, the integral involving p^(k,0) must be 
regarded as a contour integral, and the path of integration 
must be carefully chosen so that appropriate boundary 
(radiation) conditions are satisfied. 

Referring to equation 112, p. 53, shows that an inverse 

o ^ 

transform of the function Pj^(k,0) must be performed to ob- 
tain the pressure on the plate surface p(x,0,t). Comparison 

O A 

of equation 73 with 74 shows that Pj^(k,0) can be written 


(k,0) 


- (SZ* 3 


/; 


^ 12 ^ / ''2 ''2 
- (1+Mc)^ 


+ T 


22 




2 2 
- (1+Mc)"^ 


+ T 


21 


118 


provided that the square root is handled properly to yield 
the form in equation 73 when the radical changes sign. 

O 

Equation 118 shows that Pj^(k,0)has two branch points in 
addition to the simple pole resulting from the vanishing 
of the denominator. Writing the radical in the form 

/ c^ - (1+Mc) ^ ^ / 1-M^ / c + (M+1)”^ /T + (M-1)”^ 


119 
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^2 ^ 

shows the branch points to be at c = — j and c = 2 when Mach 
number M is 0.5. These branch points and their corresponding 
branch cuts are shown in figure 9 along with the pole and the 
chosen path of integration. The location of the branch cuts 
was chosen so that the radical would assume the correct form 

/N 

for real values of c as demanded by equations 73 and 74. The 
path of integration was indented below the pole for the follow- 
ing reason. Assume the integration for p(x,0,t), equation 112 

A 

was to be performed by contour integration for positive x. 

A 

With the path closed in the upper half c plane the integral 

A 

around c = infinity would vanish by Jordan's Lemma and the 
residue of the pole would be included with a contribution from 
integration along the branch cut and any other singularities pres 
ent. It can be shown that the residue from this pole contri- 

A 

butes a wave that propagates in the positive x direction; this 
wave satisfies the outgoing wave radiation condition and is 

A 

appropriately included only for positive x. 

In view of the path of integration shown in figure 9, 
the integration of equation 112 must be broken up into 
Cauchy principle parts along the real axis and integration 
along the semi-circles around the poles and branch points. 

The Cauchy principle part of the integration along the real 
axis was carried out numerically using a 10-point Runge- 
Kutta for region away from the singularity (pole) and 
a rectangular rule using the value of the function 



-2 


-] 

Path of 
integration 


I 


Branch 

cut 



I 


Figure 9.- Complex transform variable 


o> 
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at the mid-point of the interval in the region including the 

pole. The limit was taken by increasing the number of inte- 

/\ 

gration steps with a fixed c interval until the value of 
the integral was unchanged by doubling the number of steps. 
As the number of steps was increased the values of the 
function for symmetrically placed steps (with respect to 
the pole) beccime closer and closer together in absolute 
value near the singularity (as the function ^ would be for 
X = +1 and X = -1) but having opposite signs and thus 
canceling out each others contribution to the integral. The 
integral on the semi-circle around the pole was evaluated 
using the half-residue in accordance with theorem IV, p. 

530 of reference 34. The integrals on the semi-circles 
around the branch points are always zero according to 
theorem III, p. 530 of reference 34. 

The numerical values presented in this paper were cal- 
culated using the equations presented in Tables I to V. A 
sample of the functions integrated as indicated in Table II 
is shown in figures 10-13. These figures show that the 
values of the integrands are small near the ends of the c 
interval shown, indicating that integration over this inter- 
val will give correct values of the integrals. 
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TABLE I : Complete Pressure Coefficient - Including Integral 

and Singularity Contributions 


A. Complete Pressure Coefficient C^ is given by: 

1 2 j ^ 2 

1+iy-^-^l^ — § — C (x)+iTTa,Resl-yXy . ^ -R eslRes2 

2^1 (/-J^2)2| P 1 2,^4_„2^2 




where 

/s. y\ _ 

Res 1 = -« 2* ei>'^(b " 


Res 2 = ^ [(Z^)^(2=)2] (b T^ 2 +T 22 >/ ^ (b 

■ 7^ i)J 


T + 
^11 


B. To obtain: 


Cp(x) 


b, T. . 
ID 


See Table: 


III or IV 


^ "'ll ^ '^21 


R° 

Bii, 


C. The quantitiy Kq is the value of k satisfying b T^^^^ + 


0 



64 


TABLE II; The Integral Part of Pressure Coefficient C^Cx) 

A A 

A. The magnitude R(x) and phase (|i(x) of the pressure co- 
efficient (integral part) are given by 


R(x) 




<t»(x) =» tan 




C^(x) 


II-l 


II-2 


C (X) » C^(x) + i C^(x) 
P P P 


II-3 


and 


Cp(x) 


TT^ - 


1 - UB 



11 17 


TT - Q‘ 


where 


Xi (ij + i ij) 


II-4 




-20 
20 , 

c -X 


20 
20^ 


®11 27T 


(zj)^ 


<SP, 




dk, 


dk , 


II-5 


II-6 


. 7 ^ _ „ /T it cos k 

Zi - X /2 — ^ 

TT - k 


II-7 


^ /n sin k 

Zi _ X /2 -2 

TT - k 


II-8 
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TABLE II (concl.) 




3T ^ ^ ^ 

Z, cos kx + Z, sin kx I dk 


/■ *Po - r o 

J — ’’ll 


COS kx) dk 


II-9 



^Pc - r '''' c 

~ 1T ~ "^11 ^^1 ^ ^1 


6p 


r c '' 


II-IO 


B. 

To 

obtain; 

See Table 



6Pr 

X 

6 p 

^c 

' " 1 “ 

, T^j III or IV 


C. 

The 

values 

of T^j^ and T ^2 are. 



i) 

N 0 

II 

o 

^11 ^ ^12 “ ° 

II-ll 


ii) 

No 

II 

^11 = ■'ll ^12 ' ^12 

11-12 
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TABLE III; Pressure Coefficients 


6Pr «Pc 
“T"' X 


by Series 


A. Pressure Coefficients are determined by; 
2 

i) when a >0 


6p 


r ^ . '*Pc _ « ^ ^ '^12 '^22 

X X i a 


ii) when a < 0 


«Pr . , 6 ’'12 ’22 

X ^ X X b 


B. To calculate *^21 


^11 = 1 ^2 ^3 ^4 


N 


^21 == 232 233 + 4 a^ + 


n =5 

N 

Z " ®n 


n =5 


where Sq = 1 a^^ = 0 


5 M "2, 

a, 5- (1 - c ) 

^ 2X"^ 


Me 6 " 2 , 

3^ 5- (2 - c ) 

3 X^ 


^4 = 


24X 


~ x 2 w 2"2 

(1 - c2)2 _ ZMc 
4 X^ 


^n = ^n-1 -72 ^n-: 


Mc6^ Z2, _ 36^ „2: 6^ 

—5- (3 - c ) a^_3--^ M c a^_4 M d a^ 


Solution 


III-l 


III -2 


III -3 


III -4 


III -5 
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TABLE III (concl.) 


C. To calculate ’^^ 2 ' "^22* 


N 


"^12 1 + &2 + a^ + + 


^22 = ^ 2^2 


E^n 

n=5 

N 

V n a 

n 


III-6 


III-7 


n=5 


where ag = 0 a^^ = 1 a 2 = Me 


2^2 9 ''9 

, _ M^c^ _ 6^ „ _ ;2. ^ _ Mc6^ _ ;2. 

a^ ~ o o (1 c ) a. — A (2 c ) 

^ 6x^ ^ ex'^ 


III-8 


a^ = (same as for '^11' '^21^ 

2 2 

D. To determine whether a > 0 or a <0 either 

2 

i) Test a for each value of c, given M, 6, X 

ii) Use (for M < 1) 

a^ > 0 when (-1/M+1) < c<(l/l-M) 

2 

a <0 otherwise 

iii) a^ = 6^ [(1 + - c^]/X^ III-9 

= -b^ III-IO 

E. For series convergence terms are added until 

the new sum divided by the old sum is less than e. 

A typical e value is 0.001. 
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TABLE IV; 


Pressure Coefficients 
Numerical Integration 




by Runge-Kutta 


<SP- <SPc 

A. The quantities — ^ 
are as given in Table III. 


2 2 2 
a , b and the sign of a 


B. The Formulae to calculate T^^ are 


i) 








p 


'^11 

T 

12 

i ^ 





T 



dz 

1® , 

^21 

22 

hz 

0 

\ /Z=J. 



z 

U- 




- — z = 1 

- 


IV- 1 


■ N ■ 


-h. i- 


ii) 


f 1 

12 n n + 1 

Sketch of numerical integration range. 


T T 

U1 ^12 


N 


T T 

21 22 


[^N-1 ^ *-^N-2 ^ [t2][ti][to] 

iii) [t ] = ] I + §■ [A(z^) + 4A(z„ + 5-) + A(z„ + h) ] 


IV- 2 


n 


h^OhO o ®h 2®h.i 

+ -g- [A(z^ + j) A(Zj^) + + h) A(Zj^ + j) + A (z^ + j) ] 


n 2 
2 ° . h, 


+ ^ [A‘(z^ + A(z^) + A(z^ + h) A^(z^ + |)] 

4 


h"* ® 2 o h o 

+ [A(z^ + h) A^(z^ + A(z^)3 


24 


iv) 


I = 


" 1 

Lo ij 


IV- 3 


IV- 4 



N o 
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A(Zn> 


-b(Zn> 


IV- 5 




V. 


-/ ^ \ 2 “■ 

2Mc 3U ^ 1 

' , ,6,2"2 
: b(z^) ={j) c 

/I + Mcu\_ 

1 + mcu(5^)^ 

)\ 


C. Case 1; Check of previous series results 

_ o o 

Velocity profile U (z ) = z 

n n 


IV- 6 


3U 

-XT 

3z 


n 


Density profile 


PflCZn^ * constant 
3Pn 


3z 


n 


Sound speed profile 


c (z ) = 1 
n 


D. Case 2; "Blasius" profile. See Table 7.1, Reference 

O 

29 , n = 82 


_ o 

Velocity profile U(z) 


3U 

3z 


= f (n/8) 

= 8 f" (n/8) 


Density profile 


= constant 



70 



= 0 


Sound speed profile 


c(z) = 1 
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TABLE V: Derivative Entering Singularity in Complete 

Pressure Coefficient 

A. The derivative is determined from 


^ tb + TjJ = X ^ [b -f TjJ = X b 
dk dc 


dT 


11 


dc 


XTi^Sfe. X^ 

dc dc 


V-1 


where 


^ - d+MO^]"^ [(1-M^)^ - M], 

dc ^ 


V-2 


dT,, da, da- da, 

- 4- + — i + — 5- + 


N 


dc dc 


dc dc 




da 


n 

dc ' 


V-3 


and 


N 


■ A" + A 'T - A ' "T 


dc dc 


dc 


dc 


z 

n=5 


da_ 


n 


n 

— rr 

dc 


7 T 

dc 

da 


da. 


= 0 


= 0 


dc 


,r = C 

dc 


3 - <5 /o ^2, 

_ . ^ ,2 - 3= ) 


V-4 


V-5 


^®4 1 S* 


1 „2 


dc 


6 TTcd - c‘) - c^ 


i; 


= W^.vv - (n-1) (n-4) M c ^ (1 - c^) 

dc dc dc 
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w2 6^ ,, 12 , ‘^^n-3 : 6^ ^^n-4 3 '12 6^ ‘^^n-S 

Me — =• (3 - c ) — K 3M c — K- — — "M c — K- — 7i — 

dc dc dc 


I y^ 2 2 

H(H^) pn-1) ^n-2"^ ' 


20) a„_3 -3 m2 ^ - 3m" o^ ^ a^.j 


B. In the range of c considered Re {b} > 0, Im {b} = 0 


where b = 6[c^ - (1 + Me) ^3 X 
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Figure 11.- Integrand of 







Chapter IV 

Results and Discussion 
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Numerical results are presented in this Chapter for a 
plate-fluid system that can be described briefly as follows. 
Referring to figure 1/ an elastic flat thin plate lies in the 
x-y plane and is surrounded by an infinite rigid plate (baf- 
fle) also lying in the x-y plane. The elastic plate is in- 
finitely wide, has length L, is simply supported by the 
rigid plate at x = 0 and x = L, and is forced by a force 
having the distribution 

q(x,y,t) = (/Td/L ) sin ttx exp (iwt) 120. 

This force is orthogonal to all of the uncoupled plate 
modes but the first. The forcing frequency w has been 
chosen as .95 of the natural (uncoupled) frequency of this 
mode and the solution has assumed that the plate responds 
only in this mode shape. A frequency change and damping 
effect due to the fluid flowing past the plate in the up- 
per half z space are allowed for, however. (The lower 
half z space is not occupied.) The fluid properties are 
considered in two parts. The first part is a uni-direction- 
al, boundary-layer flow in the positive x direction that 
is independent of the motion of the plate. The second part 
has arbitrary direction and is dependent only on the plate 
motion. Linear equations governing the second part of the 
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fluid properties are obtained by small perturbations of the 
governing Navier- Stokes, gas-state, energy, and continuity 
equations from the first part of the flow and are solved to- 
gether with the plate equation as a coupled set of equations. 
Gradients of velocity (and temperature) , of the first part 
of the fluid flow are present near the plate associated 
with a boundary layer; these gradients enter the coupled 
plate-fluid equations as variable coefficients. The re- 
sults presented here were developed in an effort to determine 
the effects of these gradients of velocity on the response 
of the plate-fluid system. 

Pressure on plate surface 

Calculated pressure distributions on the surface of 
the plate are presented in figure 14 for three boundary 
layer thicknesses. The velocity profile was linear from 
zero velocity at the plate surface to Mach number M = 0 . 5 
at z/L = 6, and the velocity was constant at M = 0.5 for 
all larger values of z. Figure 14 presents magnitude and 
phase of pressure coefficient which is related to pressure 
by 

p(x,z,t) = PqC fL |Cp| exp (iu)t + i(j>). 121. 

The phase euigle <f> thus represents time lag relative to the 
forcing function, equation 120. Figure 14 shows that the 
pressure is largest on the surface of the vibrating plate, 

/s 

0 < X < 1, where the effect of increasing boundary layer 




M » 0.5, 6-0, 

c* - 0.493c 

0, c* - 0.995c 
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thickness 6 is a slight decrease of maximum pressure and a 
shift of the distribution toward the upstream edge of the 

A, 

plate, X = 0. Increasing 6 also increases the pressure 
both upstream (x < 0) and downstream (x > 1) of the plate. 

The phase angle <}) is nearly constant over the middle 

A. 

of the plate (x ~ .5) but varies off the plate. For large 
X the phase angle can be expressed by the linear relation 

<j) = a + b X 122. 

where a and b are constants. Inserting equation 122 into 
121 leads to 

p(x,z,t) = pQ c fL 


iCpI exp (ia) exp ji ^ ^ + x] 


123. 


Equation 123 describes a pressure wave that is moving with 
speed c' , where 

<=• = ^ ^ 24 . 

Measurement of the slope b from the curves reproduced in 
figure 14 leads to the following values for c'. 


TABLE VI: Wave speeds from pressure curves 

6 = .10 .15 .20 

1—1= 1.43 1.35 1.35 


Since b is negative these waves are moving downstream. This 
speed is somewhat greater than sound speed c' = c, but some- 
what less than the convected wave speed for this Mach niimber 
c' = 1.5 c. This result seems natural in view of the boundary 
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layer nature of the mean velocity profile. 

/N 

In the region upstream of the plate x < 0 the phase 
curve can also be represented in part by the linear relation 
equation 122. The corresponding slope of the phase curve 
is negative, as it is in the region downstream, indicating 
that "waves" are coming toward the plate from an upstream re- 
gion where no disturbances are supposed to exist. This 
result was not anticipated, and is not intuitively satisfying, 
since these waves "should" be traveling upstream, and re- 
quires further investigation. Great care was taken in the 
derivation, solution, and numerical analysis of the equations 
to avoid numerical problems , and several checks will be pre- 
sented that suggest that these unexpected waves do not arise 
from errors in the solution technique or execution. Unusual 
wave behavior has been observed previously in problems of 

/ *5 C *3 \ 

sound propagarion in ducts , ' where the anomaly was 

resolved by observing that energy is the true quantity 
whose flow must be controlled by a radiation condition, and 
that energy and pressure waves do not always propagate in the 
same direction. Investigations along these lines (among 
others) should be undertaken with regard to the waves in the 
upstream region shown in figure 14. 

Check cases M = 0 and 6=0 

To gain confidence in the methods used to obtain the 
results shown in figure 14 the limiting cases of zero 
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boundary layer thickness (6 = 0, M = .5) and zero flow 
speed (M = 0) were calculated using those methods. The 
results are shown in figure 15. No independent solutions 
have been found for comparison with these results , it 

(37 \ 

has been implied that the results for M = 0 have not 
been previously obtained. The results of figure 15 appear 
reasonable, and therefore suggest the correctness of the 
methods of this paper, for the following reasons. The 
pressure and phase angle distribution for M = 0 is symmetric 

A 

about X = .5 as would be expected in still air; the pressure 
magnitude drops off continuously with distance away from 
the plate, and in particular the phase plot indicates that 
the pressure waves are propagating away from the plate in 
both upstream and downstream directions at very near the 
speed of sound, c* = c. The curves in figure 15 for M = .5 
and (5 = 0 show that the pressure magnitude has dis-symmetry 
that would be expected for the flow situation. The phase 
curve shows that pressure waves are propagating away from 
the plate in both directions at very near the expected speeds 
of .5c in the upstream direction and 1.5 c in the downstream 
direction. 

Free Wave Speeds 

The wave speeds discussed in connection with figure 14 
are the result of two parts of the solution, namely one part 
arising from integration over the whole k range and a second 
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part arising from the contribution of the residue of the pole 
on the path of integration. In the case of figure 15 where 
M = 0 or M = .5 and 6=0 there are no poles and the wave 
nature arises entirely from the integration over k. It is 
of interest to determine the nature of the waves associated 
with the singularity and their contribution to the total 
solution. 

The transfer matrix solution, equation 72, of the govern- 
ing ordinary differential equation can be re-arranged into 
the form 


T. 


11 

-1 

1 

1 ■’^12 

SZJ 

21 

b 

in { 

°n ® 

j "'^22 


125 

the 

1 unknowns are 

0 

p and 

p. 


125, 


are the ^ , b and the Fourier amplitudes Z* of the plate 
displacement. When there is no plate motion Z* = 0 and 
equation 125 has non-trivial solutions only where 


b = 0 


(117) 


This condition also is responsible for the pole on the path 
of the k integration; the pole therefore represents a free 
wave solution for waves propagating parallel to a non-moving 
plate. In this case the solution of the fluid equation can 
be written 

p(x,z,t) = h(z) exp i k' (c't + x) 
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where the wave length of the wave is, relative to the dimen- 
sional boundary layer thickness 6*, 


Y 


wave length 
6 * 


2-nX 



126. 


and the wave speed is 

where Cq is the value of c satisfying equation 117. Values 
of c' determined for various values of y are presented in 
figure 16. These waves are downstream-propagating waves as 
indicated by the negative values of c. (See figures 11 - 13.) 
Figure 16 shows that short wave lengths (small 6) propagate 
at about the sound speed c, but long waves (large y) propa- 
gate at about the convected wave speed 1.5 c. Intermediate 
wave lengths (1 < y < 10) travel at some intermediate speed. 
Wave speeds determined from the slope of the phase curve 
of the complete solution (figure 14 & 15) are shown as the 
square symbols in figure 16, and are seen to be near, but 
somewhat below, the free wave speed (phase velocity) . The 
propagation of energy in free waves is usually associated 
with the group velocity, rather than phase velocity. Group 
velocity was determined by graphically differentiating the 
phase velocity curve in figure 16 and using the formula 

c' (group) _ c ' (phase ) _ y ^ r C* (phase) ■■ 
c ~ c ' dy ^ c ^ 


c- = i 


127, 


128. 



□ Wave speeds from complete 
solution 

M * 0.5 


Phase 

velocity- 
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The resulting group velocities are shown in figure 16 to be 
less than the phase velocity but still between sound speed 
c and 1.5 c. The sign of the group velocity is the same as 
the sign of phase velocity, indicating that for these waves 
the energy is propagated in the same direction as the pres- 
sure wave . 

The contribution of the residue (representing a free 
wave) to the pressure distribution is shown in figure 17, 
where pressure distributions calculated using only the in- 

A, 

tegral over k (the dot-dash curves) are plotted along with 
distributions including the residue term along with the 
integral. Figure 17 parts a), b) , and c) shows that 
inclusion of the residue decreases the pressure magnitude 

A 

|Cp| in the upstream region x < 0, and increases the mag- 
nitude over the rest of the x region. The changes are not 

A 

large especially on the plate, 0 < x < 1. The change of 
phase is small for 6 = .2 and 6 = .1. The difference in 
phase for 6 = . 15 figure 17b appears to be caused by a 

A 

phase shift of 360° occuring at a value of x where the pres- 
sure is small and the data points widely separated; addition- 
al calculations might remove the differences shown. For the 
smallest boundary layer value, 6 = .01, shown in figure 17d, 
the residue has a large contribution to the total solution. 
The residue contribution appears to be of an undesirable na- 
ture because it adds a downstream-moving wave in the upstream 
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region and because the agreement of this data (figure 17d) 
with the limiting case of 6 = 0 (figure 15) is better when 
the residue is not included. These results for 6 = .01 
may be due to numerical inaccuracies in the integration 
over k' since it was observed that the functions integrated 
(equivalent to those shown in figures 12 and 13, but for 

A 

6=.01) had not become negligible near the ends of the k inter- 

A 

val. The k interval was not extended because the 

series solution used was divergent outside the interval used. 
Plate response 

The response of the plate as shown be equation 112 is 
harmonic in time and has a half-sine shape over the plate 
length. Additional information can be obtained by examina- 
tion of the coefficient a^^ appearing in equation 112. The 
process of analysis leading to equation 111 for a^, when 
applied to the complete matrix equation 53, leads to the 
equation 

|[(X^ - (1 - + i B°^]j a^^ = 

129. 

The equation for a plate vibrating with an additional mass, 
m, added to its own mass Ppand with a distributed viscous 
damping is derived in appendix B, and is 


0) 


a = Q . 
m m 


(1 + ^)]| 

Comparison of equations 129 and 130 suggests the analogy 


130. 
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that the effects of the fluid on the plate can be interpreted 
as the addition of a surface mass m/p, where 



and an added viscous damping r,, where 

^ "" ^11 


131. 


132. 


Added mass and damping equivalents from equations 131 and 132 
are shown in figures 18 and 19. Figure 18 shows that the 
equivalent added mass is positive for small boundary layer 
thickness 6 and for the limit case of M = 0 but has various 
values for intermediate values .1 < 6 < 1. These results 
imply that the resonant frequency of the plate as affected 
by the added mass m can be either increased or decreased by 
the fluid effects. Similar results have been observed ex- 
perimentally. Figure 19 shows that the damping is pos- 

itive for the values of 6 used, and is about one-half per- 
cent of critical. 

Results for Blasius Velocity Profile 

The results presented so far have used the linear 
velocity profile with the series solution of the ordinary 
differential equation. The formulae are given in Table III. 
The effect of using a different velocity profile will now be 
considered. Use of a different velocity profile than the 
simple linear profile requires use of the mamerical 
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integration method of Runge-Kutta. The formulae for this 
method are given in Table IV. In order to check the numer- 

6p 

ical processes the basic pressure integrands 6p /X and — ^ 
were calculated using both series and Runge-Kutta methods 
for the linear velocity profile. A (non- typical) example 
result is shown in figure 20. In order to determine whether 
the differences shown in figure 20 are significant or not 
it would be necessary to carry through both solutions to 
pressure coefficient distributions; such a calculation has 
not yet been carried out. The differences shown in figure 
20 were the largest observed between the two methods; for 
other values of 6 the curves were close enough that they 
could not be distinguished from each other on the scale used 
in figure 20. 

Results for a Blasius boundary layer velocity profile 
are compared with results for the linear profile in figure 
21. Figure 21 shows that curves for the same value of 6/X 
for both profiles, say 6/X = 2, do not agree. However, the 
curve for 6/ X = .8 for the linear profile falls directly on 
the curve for 6/X = 2 for the Blasius profile (to the scale 
of figure 21). The definition of 6 - boundary layer thick- 
ness - is arbitrary however and is different for each velo- 
city profile. The actual velocity profiles for 6/X = .8 - 
linear profile - and 6/X = 2 - Blasius profile - are compared 
in figure 22. Considering the agreement of velocities shown 
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c = kX 


- Series solution 

- Runge-Kutta sol. 

M = 0.5 
6/X = 2.0 



c = kX 


Figure 20.- Comparison of basic integrand functions obtained 
from series and numerical solutions. 
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in figure 22, the agreement of the pressure parameters for 
the two different values of 6 shown in figure 21 is not sur- 
prising. 

The effect of the velocity profile enters the solution 
for plate response and fluid pressure only through the para- 
meters —j— and -j— . The results shown in figure 21 imply 
that a linear velocity profile can be used with an appropri- 
ate equivalent boundary layer thickness 6 to produce results 
that are the saune as if a Blasius profile had been used. 

Since the Blasius profile is known to be hydrodynamically 
stable, it can be inferred that the results presented herein 
for the linear velocity profile are valid for some hydrody- 
namically stable velocity profile. 

Constant Phase Curves 

In view of the unexpected results for wave propagation 
shown in figure 14, naimely that pressure waves appear to be 
originating in a region that should be free of disturbances, 
it is natural to seek understanding by study of the propa- 
gation of the waves. It is emphasized that the propagation 
of energy is the property that must be studied for correct 
understanding of the fluid behavior and to determine whether 
any radiation conditions have been violated. The calculation 
of energy propagation is a very complex task in a medium 
having velocity gradients and has not been carried out for 
the plate-fluid system of this paper. In particular the 



5 


2 . 


Transform variable, c = kX 


Figure 21.- Comparison of basic integrand functions for lin 
ear and Blasius velocity profiles. 
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directions of energy propagation cannot be assumed to co- 
incide with the direction of propagation of pressure waves. 

The direction of propagation of pressure waves can be 
determined using the phase plots in figures such as figure 
14. Recalling the relation between the phase curve and the 
associated wave propagation direction developed in connection 
with equations 122-124 leads to the conclusion that the waves 
flow "downhill" toward smaller values of (j>. A phase minimum 

A 

such as occurs at6=.l x=-.3 in figure 14 thus takes 
on the character of a sink for waves, whereas a maximum of 
phase has the character of a source. For non- zero boundary 
layer values the maxima and minima of phase shown in figure 
14 all occur off of the vibrating plate, whereas the maxima 
for the limiting cases of M = .5 6=0, and M = 0 all 

occur on the plate. 

Phase curves for pressure on the plate and at the edge 
of the boundary layer (6 = .2) are shown in figure 23, to- 
gether with a sketch of possible constant phase curves in- 
ferred from those phase curves. Assuming waves flow from 
large phase values to smaller phase values, figure 23 also 
shows pressure waves originating in the region upstream of 
the plate and near the plate surface, where no disturbances 
are supposed to exist. The residue contribution has not been 
included at the edge of the boundary layer, however the small 
differences shown in figure 17a on the plate surface suggest 



that the residue contribution might not be large. 







Chapter V 
Concluding Remarks 


Resume 

This thesis has described a theoretical study of the 
harmonic forced response of a thin, elastic, two-dimensional 
flat plate that is dynamically coupled to a passing air 
stream having a boundary- layer type gradient of velocity. 

The solution used a modal series solution of the plate equa- 
tion, a Fourier transform of the fluid equation in the di- 
rection parallel to the plate surface, and either Frobenius 
series or numerical integration solution of the resulting 
ordinary differential equation with variable coefficients. 
Numerical inversion of the Fourier transform leads to solu- 
tions for distributions of pressure over the plate surface 
and plate response properties. Numerical results for exam- 
ple values of the parameters are presented to illustrate 
the method of solution and the effects of the boundary 
layer on the fluid pressures and the plate response. Calcu- 
lations for the limiting cases of zero Mach number, Mach 
nximber of 0.5 with zero boundary layer thickness, and an an- 
alytical solution with zero boundary layer on an infinite 
plate that was previously solved by Miles were carried out 
to lend confidence in the methods and techniques used herein. 
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Conclusions 

The calculated surface pressure distributions showed 
that varying the boundary layer thickness changed the magni- 
tude of the pressure, however the changes were not large for 
the parameter values used. The phase of the pressure distri- 
bution showed that a pressure wave associated with the bound- 
ary layer apparently originates in the region upstream of 
the plate where no disturbances are supposed to exist. The 
explanation of the origin and nature of this wave was not 
determined. Alternate methods of solution, together with 
experimental investigation if necessary, should be used to 
explain this wave. 

Analysis of harmonic wave motion in the fluid with zero 
plate motion showed the existence of a pressure wave whose 
amplitude decreased exponentially with distances from the 
surface (and in this sense roughly analogous to Rayleigh 
waves on the surface of a solid) . At short waves lengths 
these waves travel at the sound speed c of the medium, but 
at long wave lengths they travel at the convected wave speed, 
1.5 c. 

In an effort to display the effects of the fluid on the 
plate's dynamic behavior an analogy was developed that show- 
ed the fluid effects to be equivalent to adding a distributed 
mass and damping to the plate surface. This analogy is com- 
monly used in fluid-structure dynamic interaction problems. 
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Varying the boundary layer thickness caused significant 
percentage changes in the added mass and damping, however, 
for the parameter values chosen the maximum added mass was 
only 0.1 percent of the plate mass and the added damping 
was about 0.55 percent of the critical damping. 

The effect of variations in the detail of the velocity 
profile was studied by calculating values of the fundamental 
"transformed modal pressure" parameter for a linear velocity 
profile (velocity of zero at the plate, increasing linearly 
out to the edge of the boundary layer and constant from 
there outward) and for a Blasius profile. This calculation 
showed that the same result could be obtained from each 
velocity profile provided that an appropriate value of bound- 
ary layer thickness was chosen for each profile. The appro- 
priate boundary layer thickness was shown to cause the velo- 
cities of each profile to be nearly equal in the region near 
the plate. 

Recommendations for future work 

Further topics of research suggested by these 
results include the following: 

a) Theoretical and experimental studies to clarify the nat- 
ure of the waves apparently originating upstream of the plate. 

b) Use of these methods of solution for study of panel 
flutter, as an alternative to the time-domain solution methods 
used by Dowell. 
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c) Use of these methods to search for appropriate fluid- 
panel parameters for maximizing damping of the system with 
the aim of reducing turbulence in the boundary layer and 
thus reducing drag, or with the aim of reducing noise trans- 
mitted from a turbulent boundary layer into the interior of 
an aircraft. 

d) Development of approximate procedures, based on the 
methods developed herein, that could be used in engineering 
studies instead of these lengthy calculations. 
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APPENDIX A.- Perturbation Pressure-Density Relation in a 
General Mean Flow 

In this Appendix an equation is derived relating the 
perturbation pressure to the perturbation density in an 
arbitrary mean flow. The appropriate governing equations 
are : 

(~ + u* *7) P* + p*7*u* = 0 A-1 

p* = Rp*T* A-2 

p*C (4^ + u**V)T* = -p*V‘u* + (heat) + 

V tj t 

(viscosity) A-3 

These equations are the continuity, state, and energy equa- 
tions governing the total flow quantities denoted by the 
"star". In the energy equation A-3 the terms involving heat 
flow and viscosity effects are represented simply by (heat) 
and (viscosity) because these terms will soon be dropped 
from the equations governing the perturbation quantities. 

Equation A-1 and A-3 can be combined to eliminate 
the velocity term V*u* to obtain 

(p*)^C^ D(T*) = p*Dp* + (heat) + (viscosity) a-4 

where 

D( ) = (|^ + u*-V) ( ) 

Next, using T* from equation A-2 in equation A-4, carrying 
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out the differentiations, and simplifying leads to; 

p* D(p*) = Y P* D(p*) + (heat) + (viscosity) A-5 

Note that if the effects of heat flow and viscosity are 
neglected for the total flow then equation A-5 implies that 

p* = c^p* 

2 

where c = yRT* . 

This result is ordinarily used in classical inviscid aero- 
dynamics. To obtain the relation between perturbation pres- 
sure and density, the relations 

p* = pQ + p, 

P* = Pq + 

and u* = U + u 

are substituted into equation A-5, and the resulting terms 
separated into three groups. One group of terms contains 
only basic flow quantities U, p^ , p^ , and heat flow and 
viscosity terms; these terms vanish under the assumption 
that the basic flow satisfies equation A-5. The second group 
of terms contains squares and products of the perturbation 
quantities u, p, p; these terms are dropped under the assump- 
tion that the perturbation motion is small enough that higher 
order terms are small compared to the first order terms. 

The third group of terms are linear in the perturbation 
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quantitiss and provida the desired equation i 

d(P) + 5 • 7p^ + p d(p„)/o^ = YRTp[d(p) + 5-VP(| + pd(pj)/p|,) 


where d ( ) 




• V)( ). 


A- 6 
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APPENDIX B.- Plate Equations With Damping And Added Mass 


In this Apperidix the equations governing the response 
of a plate having a mass and a damping force added to its 
surface are derived. The procedures and notations used in 
the body of this paper are used so that the equations ob- 
tained here can be compared directly with the equations 
obtained in the body of the paper for the acoustically coup- 
led plate. These equations are derived to aid in the inter- 
pretation of the acoustic-coupled equations. 

When a mass m and damping c are added to the forces 
acting on an element of the plate, the governing equation is 

2 

DV^w + (pp + m) + c 1^ = q (x,y,t) B-1 


Assuming that the plate motion is independent of y and non- 
dimensionalizing changes equation B-1 to 


4 

3 w . "P 

— — + ~rT 


8x' 


p L 2 '' 4 4 

n + + CLi ^ q(x,t)L^ 

Op 3t= “ 


B-2 


Next use the orthogonal eigenfunctions of 


9^w 


m 


3x^ 


= X w 
m m 


and the eigenfunction expansion 


w (x,t) 


= L 


an(t) w„(x) 


B-3 


B-4 


to reduce B-2 to 
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9^a 


9a_ 


m c m 

1-2 Pp (1+m/pp) 9t 


D 

m 


PpL'*(l+m/Pp) 


a = 


PpL(l+m/Pp) 




q (x,t)Wjj^(x)dx} 


B-5 


By analogy with a single degree of freedom system the 
coefficients of equation B-5 can be identified as 


Pp(l+m/Pp) " ^^“m 

- X D 



PpL (1+m/Pp) 

where is the natural (circular) frequency and c is the 

ratio of actual damping to the critical value of damping. 
Taking the force q(x,t) and the response a^it) to be 
harmonic, i.e. 



q(x,t) = q(x)e 




iut 


B-7 


and introducing equations B-6 and B-7 into equation B-5 leads 
to 


m 


B-8 


Equation B-8 is to be compared with acoustic coupled plate 


results . 



